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Abstract 

We derive a simple integral representation for the distribution of the maximum of Brownian 
motion minus a parabola, which can be used for computing the density and moments of the 
distribution, both for one-sided and two-sided Brownian motion. 



1 Introduction 

It is the purpose of this note to show how one can easily obtain information on properties of the 



distribution of the maximum of Brownian motion minus a parabola from Groeneboom (1989) In 
fact, Corollary 3.1 in that paper gives the joint distribution of both the maximum and the location 
of the maximum. In the latter paper most attention is on the distribution of the location of the 
maximum, which is derived from this corollary. The reason for the emphasis on the distribution 
of the location of the maximum is that this distribution very often occurs as limit distribution in 
the context of isotonic regression; one could say that it is a kind of "normal distribution" in that 
context. But one can of course also derive the distribution of the maximum itself from this corollary 
and at the same time deduce numerical information, as will be shown below. 

Numerical information on the density, quantiles and moments of the location of the maximum 



is given in Groeneboom and Wellner (2001), which in turn relies on section 4 of Groeneboom 



(1985) 



2 Representations of the distribution of the maximum 

Let F c be the distribution function of the maximum of W(t)-ct 2 , t > 0, where W is one-sided Brow- 



nian motion (in standard scale and without drift). Then, according to Theorem 3.1 of Groeneboom 



(1989), F c has the representation 

F c (x) = tl> x , e (p), (2.1) 
where the function ij) xc : M —> M+ has Fourier transform 



oo 



a, ,. ,., .,_ _ ir{Ai(i{)Bi(i{ + z) - Bi(iJ)Ai(i{ + z)\ 
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and where £ = (2c 2 ) 1 ^ 3 /, and z = (4c) 1 / 3 x. It follows that the corresponding density f c has the 
representation 

fc(x) = ^c(0)= f ^ c (0), (2.3) 

where (f) xc has Fourier transform 

7 m f°° us, M . ^c) 1 / 3 {MjQBi^ + g) - Bi(^)Ai'(^ + z)} 

^c(l) = J_J 4>.,c(s)ds = (2 C 2)l/3 A i(iO ' M 

The (symmetric) density g c of the maximum M of VF(i) - ci 2 , f G 1, where W is two-sided 
Brownian motion, originating from zero, therefore has the representation 

/'OO /'OO 

g c {x) = 2f c {x)F c {x) = 2(f> x , c (0)ip x , c (0) = ^ J 4>x,c( u ) dit y 0*,c(«) **, » > 0, (2.5) 

since the distribution function of M is the maximum of the two maxima one gets to the right and 
to the left of zero. Note that these two maxima are independent, since two-sided Brownian motion 
is started independently to the right and to the left, starting at zero. It is also obvious that these 
two maxima have the same distribution. Interestingly, the situation is more complicated for the 
location of the maximum! 

Note that this gives the complete characterization of the distribution of the maximum of Brow- 
nian motion with parabolic drift. The purpose of this note, however, is to show how one can deduce 
useful numerical information from this. 

The two fundamental solutions of the Airy differential equation are Ai and Bi which are un- 
bounded on different regions of the complex plane. For the purpose of computing moments, etc., 
it is easier to only work with the solution Ai, so we want to get rid of Bi. To this end we simply 
use Cauchy's formula. 

We have the following lemma. 

Lemma 2.1 Let N c be defined by 

N c = max-fPFYx) — cx 2 } , 

x>0 ' 

So N c is the maximum for the one-sided case. Then the distribution function F c = Fjsr c of N c is 
given by: 

f°° f ■,, f°° Ai (e~ i7r / 6 u) M(iu + (Acf^x) ) 
F e (x) = l- Ai (u )du-2Re{ e" l7r/6 / — i — du }, x > 0. 

One can use this representation to compute the distribution function F c in one line in, for 
example, Mathematica, and the result of this computation is shown below in Figure [TJ where we 
take c = 1/2. 



Proof of Lemma 



2.1 



After the change of variables u = (2c 2 ) 1//3 £ we get for the corresponding 
distribution function F c , still taking z = (4c) 1//3 x, 

1 f {Ai(m)Bi(m + z) - Bi(m)Ai(m + z)} 
Fc(x) = 2 AiM ( } 

Using 

Bi(z) = %M{z) - 2ie wi/3 Ai(ze- 2ni/3 ), 
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Figure 1: The distribution function Fn c , for c = 1/2. 



which is 10.4.9 in |Abramowitz and Stegun (1964)] we can write: 



1 f°° {Ai(iu)Bi{iu + z) - Bi(iu)Ai(zu + z)} 
M(iu) 



2 jo 

e 



du 



l7r / 6 / Ai(e- m/6 n + ze- 2l7r/3 ) dn-e" l7r/6 / — i A / / ' du. 2.7 

7o V ' Jo Ai(m) 



By Cauchy's formula we can reduce the first integral on the right-hand side of (2.7) to: 

J Ai (u + ze" 2 ^ /3 ) du. 

Differentiation w.r.t. z yields: 



(2. 



r "-"^ / Ai' (u + ze" 2 ^ 3 ) = -e" 2 ^ 3 Ai (z, --- /3 



We similarly have: 

1 Z" {Ai(ra)Bi(iu + z) - Bi(m)Ai(m + z)} 



CO 



Ai(ra) 



d« 



CO 



e W6 y Ai ^ e W6 u + ze 2W3 j ^ _ 



Ai (e l7r / 6 u) Ai(-m + 
Ai(-m) 



du. (2.9) 



Again using Cauchy's formula we can reduce the first integral on the right-hand side of (2.9) to: 

J™ Ai (u + ze 2in/3 ^j du. 

Differentiation w.r.t. z yields: 



(2.10) 



e 2in/3 j A y | 



du = -e 2l ^ 3 Ai (ze 2 ™' 3 ) . 



So the derivative w.r.t. z of the sum of the two integrals (|2.8|) and (|2.10|) is given by 



T / 3 Ai (ze- 2 ^/ 3 ) - e 2l ^ 3 Ai (ze 2m ' 3 ) = Ai(z). 
3 



(2.11) 



For the latter relation, see (10.4.7) in Abramowitz and Stegun (1964) 
So we get: 

J Ai (u + ze~ 2iw/3 ^ du + Ai (u + ze 2i7r/3 ) du = J Ai(u)du + k, (2.12) 
for some constant k. For z = we get: 

/*oo ^oo roo 

/ Ai(u)du+ l M(u)du = 2 I Ai(u) du = k = 2/3. 
io io JO 

Hence k = 2/3 and 



Ai (u + ze- 2i7T/3 j du + j Ai (u + ze 2iT/3 j du = j Ai(u)du + 2/3 



= 1 - y Ai(it) du. 
Note that this implies: 

lim < / Ai ( -u + ze 



The second term is given by 



-i7r/6 



Ai (e- i7T ^u) Ai(iu + 



du - e iw/6 



oo Ai ( e i7r/6^ Ai(-i« + 



(2.13) 



|y Ai (u + ze- 2iw/3 ^j du + J Ai (u + ze 2iT/3 ^j du\ = 1. (2.14) 



/o Ai M ./n 

I io Ai(wt) 



Ai(-iw) 



(2.15) 



□ 



Corollary 2.1 T/ie density of N c is given by: 
f c ( x ) = (Ac) 1 / 3 | Ai((4c) 1/3 x) - 2 Re ( e 



f 00 Ai( e - m / 6 u) Ai'(iu + (4c) 1 / 3 x) \ ] 

/6 ' ^ A ., V r — L du)},X>0. 



Proof. This follows by straightforward differentiation from Lemma 2.1 



□ 



A picture of the density /jv c , for c = 1/2, is given in Figure [2j By the remarks above, we also have: 

Corollary 2.2 The density g c of the maximum M of W(t) — ct 2 , f £ R, where W is two-sided 
Brownian motion, originating from zero, if given by 

g c (x) = 2f c (x)F c (x), x > 0, 



where f c is given by Corollary 2.1 and F c by Lemma 2.1 



A picture of the density fM c , for c = 1/2 and two-sided Brownian motion, is given in Figure [31 
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Figure 2: The density /jv c , for c = 1/2. 




Figure 3: The density Jm c = 5c of the maximum for two-sided Brownian motion and c = 1/2. 

3 Concluding remarks 

The densities of the maximum and location of the maximum of Brownian motion minus a parabola 
were originally studied by solving partial differential equations. For example, if we denote the 
location of the maximum of two-sided Brownian motion minus the parabola y = t 2 by Z, then the 



density of Z is expressed in Chernoff (1964) in terms of the solution of the heat equation 





dt 



u(t, x) 



2 & 2 



for x < i 2 , under the boundary conditions 



u{t,t 2 



def 



limu(t,x) = 1, lim u(t,x) = 0, t € K. 



x\t 2 xX—oo 

If u(t, x) is the (smooth) solution of this equation, the density fz of Z is given by 

fz(t) = \u 2 (-t)u 2 {t), x G R, 
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where (as in Groeneboom (1985)) the function U2 is defined by 

d 

u 2 (t) = lim —u(t,x). 

x^t 2 OX 

The original computations of this density were indeed based on numerically solving this partial 
differential equation (as I learned from personal communications by Herman Chernoff and Willem 
van Zwet). However, it is very hard to solve this equation numerically sufficiently accurately for 



negative values of t, since we have, by (4.25) in Groeneboom (1985) 

U2(t) ~ ci exp { — ||t| 3 — c\t\} , t — > — oo, 

where c « 2.9458 ... and c\ « 2.2638 .... 

At present the situation is drastically different, since we have much more analytical information 
about the solution and, moreover, can use advanced computer algebra packages. One only needs 



one line in Mathematica to compute the density fz, since, by (3.8) in Groeneboom (1989), fz is 
given by fz{%) = ^'M^O'K - x ) = \ u 2{x)u2{—x), where: 



4>{x) 



l 



2 2 / 3 vr 



du, 



Ai(i2-V3 u ) 

and where one can even allow the boundaries — oo and oo in the numerical integration (in Mathe- 
matica). A picture of the density fz, obtained from just using this definition in Mathematica, is 
given in Figure |4j 

However, if one wants to get very precise information about the tail behavior of the density or the 
behavior close to zero, it is better to use power series expansions or asymptotic expansions, which 
are different in a neighborhood of zero from the representation for large values of the argument. 



Details on this are given in Groeneboom (1985) and Groeneboom and Wellner (2001) 



More details on the history of the subject are given in Perman and Wellner (1996) and 

In the latter manuscript also more details on the 



Janson, Louchard and Martin-Lof (2010) 



distribution of the maximum of Brownian motion minus a parabola are given. Their results seem 
to be in complete agreement with some numerical computations, based on the representations given 
in section [2j 




Figure 4: The density fz of the location of the maximum of W(t) — t 2 , t 6 M. 
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